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Abstract 

Superconductivity is studied for a fermionic system with attractive four-body interaction. Ap- 
plying a Green function approach, the gap equation is derived. From the solution, the transi- 
tion temperature is calculated. Under the condition that the respective coupling constants are 
comparable, the transition temperature of four-fermion complexes is considerably larger than the 
corresponding BCS value. 

PACS numbers: 74.20.Fg,74.62.Yb,74.90.+n 



I. INTRODUCTION 



About fifty years ago, the development of the BCS theory^ and its overwhelming success 
in physics and application was a breakthrough in modern theoretical physics. This theory 
provides basic insight into superconductivity from a quantum-mechanical point of view. 
According to the BCS mechanism, the superconducting phase is due to the condensation of 
bosonic quasi-particles, called Cooper pairs, which are created from Fermions that experience 
an effective weak attractive interaction. An essential ingredient of the BCS theory was the 
observation^ that there could be a phonon-mediated attraction between electrons in some 
metals at low temperatures. This effective force acts at longer distances than the short- 
range screened Coulomb interaction. In the simplified version of the modeF, the attractive 
interaction is constant in a narrow energy region (of the order of the Debye energy fujjD) 
around the Fermi energy and zero outside. As the extent of a Cooper pair exponentially 
increases with decreasing coupling, a large number of Cooper pair wave functions usually 
overlap with each other. In the extreme limit of vanishing attractive coupling (72, the non- 
perturbative BCS theory becomes exact. 

The attractive interaction between fermions may originate from quite different sources. 
Since the effective attraction could be due not only to phonons but also to other bosonic 
excitations, the BCS approach has been successfully applied to many quite different research 
fields. Unfortunately, the basic physics of high-temperature superconductivity^ does not fit 
into any weak-coupling schema and therefore also not into the BCS theory. The measured 
transition temperature in cuprates, which can reach 160 K, cannot be derived from 
realistic BCS calculations. 

Regardless of this disadvantage, it is not exaggerated to state that the BCS theory 
has deepened our understanding of many-body effects in the quantum world. The non- 
perturbative character of the approach, as revealed by the highly nonlinear dependence of 
the transition temperature or the gap energy on the coupling constant g2, expresses the fact 
that the superconducting phase is to be understood as a complex many-body phenomenon, 
to which an infinite number of particular particle correlations contribute. In addition, the 
resulting Cooper pairs, which carry the superconducting current, do not exist as spatially 
separated entities but are collected into large drops, in which novel collective quantum effects 
could exist in principle. However, the BCS approach focuses on pair correlations meaning 
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that the behavior of more than two highly correlated quantum objects is described by two- 
body interactions between all possible pairs. The question arises whether this first approxi- 
mation is exhaustive for the description of many-body quantum effects. In fact, many-body 
interactions have been studied in different research areas like nuclear, atomic, and condensed 
matter physics. The main reason usually is the treatment of cluster effects created by an 
ensemble of strongly correlated particles (especially self-consistent cluster approximations^) . 
In the field of nuclear physics there is growing evidence that three-body forces exist among 
the nucleons inside atomic nuclei. In condensed-matter physics, many-body effects of trions 
were studied^ and three-body forces were shown to dominate the Cauchy discrepancies in 
the second and third order elastic constants of Copper.^ Three-body contributions to the 
interatomic potential also modify many properties of the liquid and solid phases of ^He 
atoms. ^ Even a spectacular result have been reported in the field of many-body interactiort^ 
namely that a D-dimensional Hamiltonian with two- and three-body interactions has a 
unique string-like ground state, when the strength of the three-body coupling exceeds a 
critical value. In a recent study of the Bose-Einstein condensate-BCS crossover in fermionic 
systems^, the Schrodinger equation for four attractively interacting fermions was solved by 
applying Bethe ansatz techniques. The solution reveals a robust four-fermion cluster that 
is not broken by collision and does not couple to additional fermionic states. 

The mentioned papers initiated a study of the rich and interesting quantum physics that 
refers to many-body interaction. The objective of the present paper is to contribute to 
this fascinating field by treating a BCS-like instability in an electron gas with four fermion 
coupling. 

II. BASIC APPROACH 

The BCS theory is based on the pair approximation with an effective particle-particle 
coupling that is attractive and phonon mediated. Similarly, it is conceivable that also a 
cluster of three fermions is glued together by an attracting three-body coupling. However, 
due to its fermionic character, these aggregates are not expected to form a condensate. 
One needs two or four fermion complexes for their transmutation into boson-like objects, 
which could give rise to a phase transition to a superconductor. Here, we focus on an 
instability, which is due to the conglomerate of four fermions described by an effective 
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potential v. In order to concentrate on possible instabilities in the system, both pair- and 
higher-order particle interactions are not taken into account. Therefore, we start from the 
model Hamiltonian 
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Mr,t) (1) 



+ Y1 / {^^n}^l,(ri,t)^l,(r2,t)^l3(r3,tX(r4,t) 
xv{ri,r2, rg, r4)^,, {r^, t)^}^^ (rg, t)^,^ (rs, t)^,, (n, t), 

where V{r) and m denote the crystal potential and the effective mass, respectively. Fermions 
with spin s are created [annihilated] by the operators ipKrjt) [iljs{r,t)]. Furthermore, the 
abbreviation {dvi} = dridr2dr^dr^ is used. The one-particle Green function satisfies the 
equation of motion^^ 

+ ^^r. - ^(^i) - ^^(1)1^(1, 1'; U) = 5(1 - 1') (2) 



dti 2m 

-i^^ j d2d3d4:V {1,2, 3, 1)0(1,2,3, 1; 1', 2^, 3^, 4^; t/), 

with U{1) = Us{r,t) being an arbitrary auxiliary function. In this equation, the usual 
notation is used, e.g. 2^ = {rg, % + ■^2} with e 0. As Eq. is not closed, it must be 
supplemented by an equation for the four-particle Green function and so on. This hierarchy 
is alternatively deduced from the generating functional G[X, rj] 



G(l,...,n;l',...,n';t/)=^...^^...^G[A,,] 



(3) 

X=ri=0 



by taking functional derivatives with respect to the Grassmannian source fields A and rj 
before they are set equal to zero. Expressing Eq. ([2]) by a functional differential equation 

-th' f d2d3dIV(l,2,3,A)—^ ^-T ^-T ^-=--^-^-^G[X,r]], 

J ^ ^ 577(4+) 5r/(r) 577(2+) 5A(2)5A(3)5A(4)5A(1) ^ ''J' 

the coupled set of equations for the Green functions is easily obtained by calculating addi- 
tional functional derivatives of desired order. At this stage, it is recommended to introduce 
the generating functional G'c[A,77] for the connected parts of many-particle Green functions 
via the equation 

G[A,77]=exp{GJA,77]}. (5) 



For instance, for the one- and two-particle functions, the relationship 



G(l,l';f/) = a(l,l';f/) (6) 

^(1, 2; 1', 2'; U) = ^(1, 2; 1', 2'; U) + ^(1, 1'; U)G{2, 2'; U) - G{1, 2'; U)G{2, 1'; f/), (7) 

is derived, from which it is concluded that the correlated two-particle Green function does 
not enclose the Hartree-Fock contributions. It is a general property that the correlated 
Green functions vanish, when the interaction is absent so that the hierarchy is suitably 
truncated by disregarding higher-level many-particle contributions described by correlated 
Green functions. 

Inserting the correlated four-particle Green function into Eq. ([2]), we obtain an equation, 
in which numerous terms appear that represent various kinds of four-fermion scattering. It is 
in line with our approach, which looks for a BCS-like instability for four fermion complexes, 
to neglect all two- and three-particle connected Green functions. The remaining conglomer- 
ate of one-particle Green functions is collected into an effective function Go{l, 1'; U) [in the 
pair approximation, which is the basis of the BCS theory. Go is given by the Hartree-Fock 
Green function] so that we obtain 

1'; U) = Go(l, 1'; U) - ih^ j dld2dMlGo{lJ\ f/)V(T, 2, 3,4) 

xGe(T,2,3,4;l',2+,r,4+;?7). (8) 

This result is complemented by the equation of motion for the eight-point Green function, 
which is derived from Eqs. (jlj) and ([5]). The handling of this rather complicated integral 
equation is guided by the well established pair approximation^iii^ that leads to the BCS the- 
ory. In accordance with the BCS reasoning, we focus on the homogeneous Cooper channel, 
in which the connected four-particle function is self-consistently determined from 

Ge(l,2,3,4;l',2',3',4';f/) =in^y" rfTd2d3c/4Go(lJ; f/)^(T, 2, 3,4) (9) 
xG(2, 2+; f/)G(3, 3"-; f/)G'(4, 4+; U)G,{X % 3, 4; l', 2', 3', 4'; U). 

A diagrammatic representation of the scattering contribution is shown in Fig. 1. The basic 
set of Eqs. (IE]) and ([9]) extends the BCS theoryii»^ to fermions coupled via a four-body 
potential. The solution of these equations is facilitated by a Fourier transformation in the 
time domain. Introducing the Matsubara frequencies {zy = mu / {hP) + fi/h, u = ±1, ±3, . . . 
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FIG. 1: Diagrammatic representation of the Cooper channel. The wavy hnes represent the four- 
body interaction. 



and uJn = iTm/{h(3) + 4:fi/h, n = 0, ±2, . . . with fi being the chemical potential and (3 
l/ksT), the Fourier transformation is defined by 
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a(l,2,3,4;l',2',3',4';t/)= — ^ J] J] KUK}; K^j^) 

X exp [-iz^^ti - izy^t2 - iz^.^h - - z^^ - z^^ - z^^)t^ 



X 



exp [izyit[ + iz^^'t'^ + iz^'^t'^ + i(a;„ - - z^'^ - z^^%] 



(10) 



where the spin and spatial coordinates are expressed by {xi} = xi,X2,X3,X4 with xj = Vj, sj 
and {zui} = Zy^, z^^, Zi^^. When treating the Fourier transformed basic Eqs. ([8]) and Qj, the 
following functions occur 

G,{{x.y,{x'M^y^}M = (^) J]a({x.};{x:}|KJ;K.}|c..), (11) 
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a({x.}; {x'^M =(-.] J2 ^^(M; WM^uA; {z.^}\ 

A({a;.}; {x[}\u^) = v{{r,})v{{r:})G,{{x,}; {x[}M. 



(12) 



(13) 



The integral equation for the gap function A in the representation with Matsubara frequen- 
cies is straightforwardly obtained from Eq. iQ 
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A{{xi}; {x'i}\ujn) = -iffv{ri,...,r4) j {dxj} ( — J ^ G'o(a;i, x-, 2;^JG(x2, ^..J 



xG{x:i,x^,Zy.^)G{xi,xj,Un- Zy^ - Zy^ - 2:1,3) A({a;J; {a;.}|a;„). 



(14) 
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To derive the corresponding Fourier transformed version of Eq. ([8]), we account for the 
symmetry property 

{z,':]M = -G:{{x',y, {x,}\{z:,}; (15) 

which leads to an alternative formulation of Eq. f[T^ . Inserting this gap function into Eq. 
and performing a spatial Fourier transformation by using a definition similar to Eq. fllOl) . 
we obtain 

X f ^ J X] X] G{ki, Zy^)G{k2, Zy^)G{q - k - ki - k2;ujn - - Zy^ - Zy^), 
where the new gap function is given by 

A(g, ^n) = Y^Y.Y. ^".). (17) 

{s^} {fcj {fca 

In the derivation, it was considered Gss' = Sss'G, which dictates the spin dependence of the 
gap function. The basic Eq. (fT6!) for the one-particle Green function provides an obvious 
extension of the BCS singlet state description.— Similarly, we obtain for the gap equation 
the following Fourier-transformed version 

A{{ki}; {ki'}\q,uJn) = -2^^ ^f(fci - kj, ki- kj+k2- % 

, fci - fci + fca - % + fcg - %) ( ^ ) V Go{kj, z^JG{k2, Zy,,)G{k^, z^^) 

^^^^ W 

xG(q - % - % - %,t^„ - Zy^ - Zy^ - Zy.^)/\{{kj}] {ki}\q,ujn). (18) 

In the spirit of the BCS theoryiiii^, Eqs. (1161) and f|T8l) are solved in the limit q = and 
n = 0, which is an approximation concerning the effective bosonic excitation composed of 
four particles. In addition, it is assumed that the effective attraction is given by a contact 
interaction, which is confined to a narrow shell of width 2huJD around the Fermi energy. 
Introducing the new abbreviations k2 + k^ = q and Zi,^ + z^,,^ = Um, we obtain the final set 
of coupled equations 

G{k,Zy) = Go{k,Zy) + Ah^Gl{k,Z^) ( ) S2u{q,UJm)G*{k + q, Zy + UJm), (19) 
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q,m 



-ih^gi f^j f^j ^G'o(fc,^^)G*(fc + q,z^ + a;„) = 1, (20) 

^ ' q,m ^ ^ fc,i^ 

M(g,u;„) = f^j ^ G'(fci,2;,,JG'(g - - z^J, (21) 

fci,i/i 

which generahzes the BCS theory to a fermionic system with four-body interaction. These 
equations allow a detailed description of the superconducting phase of the considered cluster 
model. Formally, Eqs. (fT9|) and (120!) reduce to the basic BCS equations, when the function 
u{q,ujm) is replaced by -ihl3Sm,oSq,o- 

III. TRANSITION TEMPERATURE 

The internal structure of the four-particle cluster is mainly accounted for by the function 
u{q,u!n), which has a bosonic character. To study its influence in more detail, let us treat 
the transition temperature T^, at which the gap energy A vanishes. For the calculation of 
Tc, the self-consistent one-particle Green function G is given by Gq, for which, in analogy 
to the BCS approach, the simple expression 

«»<*-^"> = s:b^- p^' 

is adopted, denotes the kinetic energy h'^k'^/2m. For the function u{q,uJn), we immedi- 
ately obtain 

u{q,uji) = -y ; , (23) 



fc 

with the Fermi distribution function np, iuJi = 27iil/Pc, Pc = ^/ksTc, and = £k — yU- 
Inserting this result into Eq. fl20l) . we obtain the gap equation 

9a^Y1 i^F{^ki+q) + npi^ki) - 1] [nF{^k2-q) + ^^(^fca) - 1] 
fci,fc2 q 

^ nsi^ki+q + ^fci) + nB{^k2~q + ^fca) + 1 ^ ^ (^24) 

^fcl+q + Cfci + ^k-z-q + ^fc2 

in which the Bose distribution function ub and the four-body coupling appear. Let us 
compare this result of the BCS-like approach for four fermions, which are held together by 
a four-body phonon-mediated attraction, with the gap equation of the BCS theory 

^ 2n^(a) - 1 
92}_, = 1, (25) 



the solution of which is given by the well-known formula 

which applies under the conditions that there is an attractive pair interaction and that the 
inequalities /i ^ huo, /3c^d/2 ^ 1 are satisfied. A^o denotes the density of states at the 
Fermi energy. To solve the gap Eq. ( !24ll of the four-fermion cluster, it is considered that the 
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FIG. 2: Transition temperature for the g'4 approach (upper curve) compared with the BCS re- 
sult (lower curve) for hujo = 0.02 eV. The effective dimensionless coupling strengths g2NQ and 
54A^o^-d/4 are assumed to be equal. 

energies e^-^ and e^^ are confined to a narrow interval around the Fermi surface so that the 
main contribution comes from q = 0. Adopting this approximation, we arrive at 

?±Ni I dE,iE, tanh (Ml) tanh (Mi) -"^"''^'^^j + g"'^'^-' = 1. (27) 

The numerical solution of this equation for characteristic coupling strengths of Cooperons 
is shown in Fig. 2 by the upper curve and compared with the BCS result (lower curve). For 
the comparison, the effective pair and four-body coupling constants were set to be equal. 
In this case, the transition temperature of four-fermion clusters is considerably higher than 
the BCS value. Whether the predicted high-transition temperature is realistic depends on 
possible values of the four-fermion coupling strength. The study of this problem requires 
further work that goes beyond the scope of our paper. 
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IV. SUMMARY 



Guided by a Green function approach to BCS superconductivity, a superconducting phase 
transition has been identified in a fermionic system with a weakly attracting four-body in- 
teraction. The basic building blocks of the superconducting phase are quasi-particles that 
consist of four fermions hold together by an effective attraction. Special results have been 
obtained for the transition temperature Tc by solving the gap equation. For comparable cou- 
pling strengths, the transition temperature of the four-fermion cluster model is much higher 
than the BCS values. This kind of superconductivity should be observable (preferentially 
by Andreev refiection) for strong coupling strengths within the four-fermion conglomerate. 
Whether the agglomerated attraction of four fermions can reach characteristic BCS values 
is a question that will decide future studies. 
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